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Abstract 



Recently, there have been some interesting developments involving off- 
forward parton distributions of the nucleon, deeply virtual Compton scat- 
tering, and hard diffr active vector- meson production. These developments 
are triggered by the realization that the off-forward distributions contain in- 
formation about the internal spin structure of the nucleon and that diffractive 
electroproduction of vector mesons depends on these unconventional distri- 
butions. This paper gives a brief overview of the recent developments. 
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I. INTRODUCTION 



Parton distributions were first introduced by Feynman as the phenomenological quanti- 
ties describing the properties of the nucleon manifest in high-energy scattering |]]. When 
a nucleon travels near the speed of light, it can be viewed as consisting of a beam of non- 
interacting partons that are massless and collinear. Parton distributions are the densities 
of partons as a function of the fraction x of the nucleon momentum they carry. With the 
advent of Quantum Chromodynamics (QCD), Feynman's parton model, modulo logarithms, 
can be justified in a rigorous field theoretical approach 0. Meanwhile, parton distributions 
can be expressed as the matrix elements of the bilocal light-ray operators in the ground 
state of the nucleon ||. Only in the light-like gauge, A + = 0, and light-cone quantization, 
does their physical interpretation become apparent. 

When studying the physics of parton distributions themselves, one of course can no 
longer regard the nucleon as consisting of free partons. Calculating these distributions 
requires a detailed knowledge of the nucleon wave function in quark and gluon degrees of 
freedom, and undoubtedly is one of the most challenging problems in nonperturbative QCD. 
At present, the only fundamental approach we know is the lattice QCD method M. On the 
phenomenological side, however, the unpolarized parton distributions have been extracted 
with good precision from various high-energy scattering data taken in the last thirty years 
0-0. The empirical parton distributions not only provide the necessary input for calculating 
high-energy scattering cross sections, e.g., of top-quark production, but also contain valuable 
information about how the nucleon is made of its constituents. 

Recently, there has been much theoretical activity surrounding the so-called off-forward 
parton (or non- forward, off-diagonal) distributions (OFPD's). The new distributions gen- 
eralize those introduced by Feynman: They characterize certain properties of the nucleon 
exhibited in a class of high-energy scattering; they reflect the low-energy internal structure 
of the particle. However, there are also important differences between the off-forward and 
forward distributions. The OFPD's in general cannot be regarded as particle densities, 
but rather their physical interpretation is given in terms of probability amplitude. As we 
shall see, they do have simple physical significance in light-cone coordinates (or the infinite 
momentum frame). 

The most natural appearance of the OFPD's is in non-forward high-energy processes. 
As far as we know, such processes were first considered by Watanabe || and by Bartels 
and M. Loewe ||. In Ref. |§, a general virtual-photon Compton scattering off a nucleon 
target was studied in operator product expansion (OPE), an integral representation of the 
amplitude was obtained, and the off-forward matrix elements of twist-two operators were 
defined. The study, however, mentioned little about physical motivations and the practicality 
of the process; in particular, significance of the integral representation was obscure. In Ref. 
||, diffractive production of Z° and photons was studied in perturbative QCD and the 
properties of the non-diagonal gluon ladder were investigated. However, no notion of an 
off-forward distribution of the nucleon target was apparent in the paper. Following M, 



a set of off-forward evolution kernels were presented by Gribov, Levin, and Ryskin |T0 
These kernels were considered as those for the scattering amplitude. The role of off-forward 
distributions in hard diffractive processes had remained obscured in relatively recent studies. 



In Ref. [11 1, diffractive electroproduction of J/tp was analyzed by Ryskin in the small x 
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region. The production amplitude was related to the usual Feynman gluon distribution. 
The off-forwardness was taken into account by a phenomenological form factor. In Refs. 
r^T4[, Ryskin's analysis was generalized to light vector-meson production in which the 



meson structure is accounted for by the leading-twist light-cone wave function. 

To the author's knowledge, one of first examples of OFPD's was introduced by Dittes 
et al. fl5 |, and was called "interpolating function." Ironically, the main motivation for 



introducing such an object was not from studying any off-forward process; rather, it was 
considered as an example of the physical quantities whose evolution kernel interpolates be- 
tween the two well-known cases: the DGLAP and ERBL evolutions JTB|-rS|. A serious study 
of the OFPD in a physical process was made by D. Miiller et al. |19j in which they showed 
that it contributes to the general two virtual-photon Compton scattering — the same process 
considered earlier by Watanabe. They also derived the evolution equations for the OFPD. 
Note that in Ref. a general off-forward amplitude (or single-particle Green's function) 
was introduced in studying hadron helicity flip in hard scattering. However, because the 
correlation of two quark fields in the object is off the light cone, it does not have a simple 
parton interpretation as the OFPD's do. 

The recent wave of interest in the OFPD's was generated from two new developments. 
The first development is the recognition of the importance of the OFPD's in describing the 
internal structure of the nucleon pi|. I* 1 studying the spin structure, it was found that the 



fractions of the nucleon spin carried by quarks and gluons can be determined by the form 
factors of the corresponding energy-momentum tensors. These form factors can be obtained 
from the moments of the OFPD's which were independently introduced in In the same 



paper, deeply virtual Compton scattering (DVCS) was proposed as a practical means to 
measure the OFPD's. In DVCS, a highly virtual photon scatters off a quark in a nucleon 
target, which then emits a real photon and returns to form a recoil nucleon. Although 
the two-photon process studied earlier by Watanabe and Miiller et al. contains DVCS as 
its special kinematic limit, the latter has important experimental advantages and a special 
theoretical status. Indeed, DVCS is more delicate to treat theoretically because of an extra 
light-like vector among the external momenta. 

The second development is the realization that hard diffractive electroproduction of 
vector mesons requires the use of OFPD's p2|j2l[| . The first serious studies in this direction 
were made by Radyushkin |23[, Hoodbhoy f2H}| , and Collins et al. [p5] . In Ref. p3| , an off- 
forward gluon distribution - most important in the low x region - was introduced and used 
to calculate vector meson production in deep-inelastic electron scattering. The result is free 
of the assumptions made in Ref. ]T2|] about the dominance of the absorptive contribution 
and the use of the forward gluon distribution. Hoodbhoy fM\ used the off-forward gluon 
distribution to calculate heavy quarkonium production investigated earlier by Ryskin pT 



He obtained a similar result as Radyushkin's. In a comprehensive paper, Collins, Frankfurt, 
and Strikman showed that a factorization theorem exists for general hard diffractive meson 
production in deep-inelastic scattering 
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The factorization theorem explicitly involves 
the OFPD's. 

Recently, many works have been produced in studying off-forward parton distributions, 
deeply virtual Compton scattering, and diffractive vector meson production. It is the goal 
of this paper to review these interesting developments. In section II, off-forward distribu- 
tions are introduced from different perspectives and various definitions in the literature are 
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compared. In section III, we make a digression to discuss the spin structure of the nucleon. 
The role of the OFPD's in this interesting subject is exposed. In section IV, we summarize 
our present knowledge about the OFPD's at hadron mass scales and their evolution to high- 
energy scales. In section V, we consider processes in which the OFPD's may actually be 
measured: deeply virtual Compton scattering and hard diffractive vector-meson production. 
Section VI contains the summary and outlook. 



II. DEFINITIONS AND INTERPRETATIONS 

Since the off-forward distributions contain rich information about the nucleon structure, 
it is instructive to examine their definitions and interpretations from different angles. In 
this section, we will first introduce the off-forward distributions from the form factors of 
twist-two operators and then we explore their significance in light-cone coordinates. Finally 
we compare the OFPD's with other equivalent definitions in the literature. 



A. Elastic Form Factors of Twist-Two Operators 

From the point view of the low-energy nucleon structure, it is perhaps most interesting 
to consider the off-forward distributions as the generating functions for the form factors of 
the so-called twist-two operators. Recall that the matrix elements of the electromagnetic 
current in the equal momentum states are entirely determined by symmetry, whereas those 
in the unequal momentum states define the (Dirac and Pauli) form factors which contain 
such interesting information as charge radius and magnetic moment of the nucleon. The 
following tower of twist-two operators represents a generalization of the electromagnetic 
current 

of-"" = v 111 iT" 1 7 Mn V 9 , (i) 

where all indices are symmetrized and traceless (indicated by (■■■)) an d V= {V — V)/2. 
Technically, these operators transform like (n/2,n/2) under Lorentz transformations. Our 
interest in such operators is not accidental. They appear, for instance, in the operator 



production expansion of the two electromangetic currents in deep-inelastic scattering p6 
Thus although these generalized currents do not couple directly to any known fundamental 
interactions, they can nonetheless be studied indirectly in hard scattering processes. 

Since the operators for n > 1 are not related to any symmetry in the QCD lagrangian, 
their matrix elements between the equal momentum states, 

(P\0^-^\P) = 2a n (/i)P (An • ■ -P^ ] , (2) 

contain valuable dynamical information about the internal structure of the nucleon. In this 
paper, we use the covariant normalization (P\P) = 2E(2tt) 3 5 3 (0) for the nucleon state. The 
fi dependence of the above matrix elements signifies the dependence on renormalization 



scale and scheme ||26fl . The quark distribution q(x, fi) introduced by Feynman has a simple 



connection to the above matrix elements: 
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1 

dxx n ~ q(x, /i) = a n (fj) , (3) 

where q(x, fi) is chosen to have support in (—1,1). For x > 0, q(x,fi) is the just density 
of quarks which carry the x fraction of the parent nucleon momentum. The density of 
antiquarks is customarily denoted as q(x,fi), which in the above notation is —q(—x,fj,) 0. 

Just like the form factors of the electromagnetic current, additional information about 
nucleon structure can be found in the form factors of the twist-two operators when the matrix 
elements are taken between the states of unequal momenta. Using Lorentz symmetry and 
parity and time reversal invariance, one can write down all possible form factors of the spin-n 
operator 

(P'\0^-^\P) = U(P')j^U(P) E A qnt2i (t)A^ ■ . . A^+ 1 P At2l+2 ■ --P^ 

j- n — 1 j 

+ U{P') a{ ^ a U{P) E B qn ^{t)^...^^p^ +2 ...p^ ] 

+ C qn {t) Mod(n + 1, 2) -Lu(P>)U(P)A^ ■ ■ • A"»> , (4) 

where U(P') and U(P) are the Dirac spinors and Modin + 1,2) is 1 when n even, when n 
odd. The four-momentum transfer is denoted by A = P' — P and its invariant by t = A 2 . 
The average nucleon momentum P = (P' + P)/2 is also a useful variable. For n > 1, even 
or odd, there are n + 1 form factors. C qn (t) is present only when n is even. 

As the forward matrix elements a n , the above form factors can be used to define a new 
type of parton distributions — off-forward distributions. To accomplish this, we introduce 
a light-like vector n M , which is conjugate to P in the sense that P ■ n = 1. Write P = 
p+ (M /2)n, where M = M 2 — t/4 and p is another light-like vector. Contracting both 
sides of Eq. with n Ml ■ • • n Mn , we have 

n w ■ ■■n, n (P'\0^\P) = U(P') iU{P)H qn {^t) + U(P') ^ a U(P)E qn (£,t) , (5) 



where 

H qn (U) = E ^n,2,(t)(-20 2j + Mod(n + l,2) C qn (t)(-20 n , 

i=0 

E qn {^)= E ^ n , 2j (t)(-20 2l -Mod(n + l,2) C qn {t){-2£) n . (6) 

i=0 

Here we have defined a new variable £ = —n- A/2. Clearly, the £ dependence of H qn and E qn 
helps to distinguish among the different form factors of the same operator. Now one can 
introduce the off-forward parton distributions (OFPD's) H(x,^,t) and E(x,£,t) through 
their moments 

f_ dxx n - l E q (x, £, t) = E qn {^ t) , 

f 1 i dxx n - 1 H q (x,Z,t)=H qn (U) ■ (7) 
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Since all form factors are real, the new distributions are consequently real. Moreover, be- 
cause of time-reversal and hermiticity, they are even functions of £. For simplicity, we have 
not shown here the renormalization scale dependence. Frequently we will also omit the t 
dependence. 

The new distributions are more complicated than the Feynman parton distributions 
because of their dependence on the momentum transfer A. As such, the OFPD's contain 
two more scalar variables besides the Feynman variable x. The variable t is the usual t- 
channel invariant which is always present in a form factor. The £ variable is a natural 
product of marrying the concepts of the Feynman distribution and form factor: The former 
requires the presence of a prefered momentum along which the partons are predominantly 
moving and the latter requires a four- momentum transfer A; £ is just a scalar product of 
the two momenta. Some insights about the dependence of the OFPD's on these two extra 
variables will be offered later in Section IV.A. 

In QCD, there are five additional towers of twist-two operators besides that in Eq. ([!]): 



q 


— <->(Mi 

= t/> q iV ■■ 


<-+Mn-l 

■iV 




U qT 


— 

= Vv v ■ • 


■i V 




a 


= F^ ia i V 


■■■iV 


A a 


3 


= F { ^ a i V 


■■■iV 




/npi—finaP 


= F iflia i V 


■■■iV 


_1 



The corresponding off-forward distributions can be labelled by (Jl q (x,£), E q (x,£)), 

(H Tq {x,0,E Tq (x,£)), (H g {x,0,E g {x,0), (H g (x,£),E g (x,£)),an.d (H Tg (x, £), E Tg (x, £)), 
respectively. 



B. OFPD's in Light-Cone Coordinates and Gauge 

The physical significance of parton distributions in high-energy processes is apparent 
only in light-cone coordinates and light-cone gauge. To see this, we sum up all the local 
twist-two operators into a light-cone bilocal operator and express the parton distributions 
in terms of the latter, 



iXx 



P' 



P 



1«, 



1 irr^n A 

1 — , ;V>. " ^ V U{P) 



(9) 



The light-cone bilocal operator (or light-ray operator) arises frequently in hard scattering 
processes in which partons propagate along the light-cone. In fact, the Taylor-expansion of 
this operator along the light-cone leads us immediately to the twist-two operators considered 
in the last subsection. The parton distributions are most naturally defined in terms of the 
matrix elements of the bilocal operator. In this context, the Feynman x is just the Fourier 



6 



variable of the light-cone distance. The Lorentz structures in the second line in the above 
equation are independent and complete. 

In ordinary coodinates, the correlation of two quark fields in the light-ray operator is in 
both space and time. As such, the physical interpretation is quite complicated. However, 
in light-cone coordinates, the physical content of the bilocal operator becomes transparent. 
To see this, let us recall that the light-cone coordinates are defined as, 



x 



± 



-^(x°±x 3 ) 



1 

•jC _|_ ii^pC - t£ ^ • 



(10) 



In this new system of coordinates, the role of time is played by x + . In the same fashion, 
introduce the light-cone Dirac matrices ^ and the projection operators, P± = |7 =F 7 ± . A 
Dirac field if) can be decomposed into a sum of if> + = P + if>, the independent degrees of 
freedom ("good"), and if)_ = P-ip, the dependent degrees of freedom ("bad"). In light-cone 
quantization, the independent components of the Dirac field obey the following commutation 
relation f2"7j 



if) +a (x),ip+p(y) 



\x~ r =y 



1 

V2 



P +a p5(x -y )S 2 (x ± -y±) 



(11) 



where [ ] + denotes anticommutation. The above relation can be solved by a Fourier expansion 

dk + d 2 k\ ^ /. , r» * i( x -k+-x ± -k ± ) 



+ d{(k + ,k x )v x (k)e^- k+ - s ^) , (12) 

and similarly for if> + . The quark (antiquark) creation and annihilation operators, b\ k (d\ k ) 
and b\k (d\k), obey the commutation relation 



bx(k + , k ± ), b\,{k' + , k> ± )} , = 2k + 5 xy (27rfS(k + - k' + )5 2 {k L - k' ± ) , 



dx(k + , k ± ), d[,(k' + , k\) = 2k + 5 X y(27i) 3 5(k + - k' + )5 2 (k ± - k' ± ) . 



(13) 



In the presence of a gauge field, the above is unchanged if the light-cone gauge, A + = 0, is 
chosen. The same gauge choice eliminates the gauge link in the light-ray operator. 
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FIG. 1. Diagrams showing £ symmetry of the off-forward parton distributions: a) x > £; b) 

£ > x > c) x < 

Now we can look at the physical content of the off-forward distributions in the light- 
cone coordinates and gauge. Since they are even in £, we first restrict ourselves to £ > 0. 
Substituting the light-cone Fock expansion into Eq. (|9p, we have 



F q (x,0 



1 



d 2 k\ 



E 



in) 3 x 



2p + V J 2 ^\x 2 -f 2 |(2 

'P' b\ ((x - 0P + , k ± + A ± ) b x ((x + 0P + , k 



x { IP' 



Pj, for x > £ 

((-x + £)P + , -*± - A ± J 6_ A ((x + £)p+, £ ± ) | P) , for £ > x > -£ (14) 
4 ((-^ - 0/> + , fei + A ± ) d x ((-x + £)P + , P) , for x < -£ 



P 



where V is a volume factor. The distribution has different physical interpretations in the 
three different regions. In the region x > £, it is the amplitude for taking a quark of 
momentum k out of the nucleon, changing its momentum to k + A, and inserting it back 



S 



to form a recoiled nucleoli. In the region £ > x > — £, it is the amplitude for taking out a 
quark and antiquark pair with momentum —A. Finally, in the region x < — £, we have the 
same situation as in the first, except the quark is replaced by an antiquark. The first and 
third regions are similar to those present in the ordinary parton distributions and hence will 
be called the DGLAP region. The middle region is similar to that in a meson amplitude 
and hence will be called the ERBL region. 

The £ symmetry is an important property of the OFPD's. To illustrate the parton 
content of this symmetry, one can work out the equivalent of Eq. fll4|) by assuming £ is 
negative. The result is best shown in Fig. 1, where the quark and antiquark creations and 
annihilations are indicated explicitly. 

From Eq. fll4|) , the off- forward distribution is very much like a quark- nucleon scattering 
amplitude. Use Auh.wh 1 to label the helicity amplitudes, where H, H' are the helicities of 
the initial and final nucleons and h, hi are those of the created and annihilated quarks. Time- 
reversal and parity invariance provide the constraints Anh,H'h' = Ah'h^hh an d Auh.H'h' = 
A-H-h-H'-h 1 - Therefore, there are total of six independent helicity amplitudes which imply 
four additional twist-two off-forward distributions 



F Tg (x,0-- 



dX , 



2tt 



P' 



X 



-n 



P 



1 _, l ™ A • n 



dX 
— < 

2tt 



P' 



A \ 

-n I n^a 



1^ 



H Tq (x,0 -U(P')n^ ± U(P) + E Tq (x,0 -U(P' 



M 



U{P) 



(15) 



Their parton interpretation is similar to that of F q (x,£), except that for F q (x,£) there is a 
sign depending on the quark helicity and for there is a quark helicity flip 



Off-forward gluon distributions can be introduced in the similar fashion. Consider first 
the gluon-nucleon helicity amplitudes. Symmetries limit the number of independent ones 
to six, and hence there are six leading-twist gluon distributions. Here we quote only the 
helicity-independent ones: 



dX 

2xJ 2tt' 



P' 




P ) 



1 1 ia^n A 

H g (x,0 -U(P') jU(P) + E g {x,Z) -U(P') 2J ^ V U{P) 



(16) 
-F G {-x,0- 



where the gauge link is not shown. Because of the Bose symmetry, Fq(x, £) 
The interested reader can find the other twist-two gluon distributions in f28|. 

In light-front coordinates, only A± are the independent degrees of freedom. They obey 
the commutation relation [ET7H 



[A i (x),A j (y))\ :B + = y+ = -<%<f(xL - y±)e(x - y ) 



(17) 

where e(a — b) = 6(a — b) — 6(b — a). The commutation relation can be solved by the following 
expansion 
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Ai(x ,x±) 



rlh + rl 2 h, 



+al(k + ,k ± )e*(\,k)e 



i(x k + —x±-kj_) 



(18) 



where a\(k) and a\(k) are the gluon creation and annihilation operators. Using this expan- 
sion, the physical content of the gluon distributions can be worked out immediately 



d 2 k\ 



E 



Ap+V J 2 J\x 2 -e 2 |(2vr) 3 a 



P' 
P' 



(x - C)P + , k± + A ± ) a x ((re + 0p + , k± 
{~x + 0P + , ~k± - A±) a_ A ((a; + 



P 



for x > £ 
P) , for £ > x > 



(19) 



In the region x > £, it is the amplitude for taking a gluon of momentum out of the nucleon, 
changing its momentum to k + A, and inserting it back to form a recoiled nucleon. In the 
region £ > x > 0, it is the amplitude for taking out a gluon pair with momentum —A. 

In an actual calculation of the distributions, the vacuum subtraction in matrix elements 
is always implied. As such, it is convenient to insert a T product between the two fields at 
light-cone separations. The difference between a T product and an ordinary product is just 
a constant 



TV+(0)V+(An) = ^i(0)V+(An) - 6(\n°) UJ.(0), V+(An) 



(20) 



The constant does not contribute to any physical matrix element. Therefore, a light-cone T 
product contains only dispersive contributions. 



C. Other Notations and Definitions 



The notion of OFPD's we discussed above are essentially that of Refs. [H5| and |^TJ. The 
variables x and £ here are the same as the variables t and r in fl9f . Other names (non- 
forward distributions, off-diagonal distributions, double distributions) and notations exist 
in the literature. 



Radyushkin introduced the non-forward distributions ^(X) in Refs. [P3ip0| . To see their 
relationship to the OFPD's, one can write Eq. (0) for an arbitrary light-like n 



iXk-n I p 

2tt \ 



da n-A(an) , ( A 

V2 U - n 



n ■ k n ■ A 



n-P 2n-P n-P 



1 u(p') iu{p) 



+ E, 



n ■ k 

'/ \ Z 75' 



n-P In- P) n-P 



„io -^n^A 

2M { } ' 



P 



(21) 



One can also shift the argument of the string operator using Heisenberg's equations of 
motion. Now if one chooses 
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n-P = 1, X = k-n, C 
the variables are related according to 

X-C/2 



-A ■ n 



(22) 



x 



X 



l-C/2 



C/2 . 

l-C/2 ' 



c 



2£ 



The matrix element of the string operator becomes 



(23) 



fo(0) J lVe- ia ^ da n ' A[an) %{\n) 



P) = F q (x,Z,t) 



—e^ x ( P 
2tt 

jrr^n A 

(l + 0H q (x,0U(P') iU{P)+ (l + 0E q (x,0U(P') V 



U{P) 



Compared with Eqs. (9.1) and (9.2) in Ref. [^0|1 , we have 

(1 + £)H q (x, = ^(X)#(0 < X < 1) - jf(c - X)0(-1 + c < x < c) 



(24) 



(25) 



The advantages of H(x, £) over JF^(X) are: 1) it treats the initial and final nucleon symmet- 
rically and thus £ symmetry is obvious; 2) it has a simple connection with local hermitian 
operators. However, the X variable in JF^(X) has the advantage of measuring one of the 
initial parton momenta according to the longitudinal momentum of the initial nucleon, anal- 
ogous to the forward case. 

Collins et al. introduced the off-diagonal distributions f(xx,x 2 ,t) p5| . They also chose 
the n vector conjugate to the initial state momentum. The relations of their variables to 
ours are, 



x 1 



X 



X 2 



x-£ 



{x 1 + x 2 )/2 



{x x - x 2 )/2 



l-(xl-x 2 )/2 , l-{xx-x 2 )/2 



(26) 



Here x\ is the same as Radyuskin's X. The off-diagonal distributions are related to the 
OFPD's by 



fg/p(x 1 ,X 2 ) = Fg{x,£). 

On the other hand, for gluons, one has 

xix 2 f g / p (xx,x 2 ,t) = 2xF g (x,£) . 



(27) 



(28) 



Again, the £ symmetry is not obvious in these off-diagonal distributions. 

Radyushkin also introduced what he called double distributions F(x,y), where x and y 
are the fractions of the initial nucleon momentum P and the momentum transfer A carried 
by the active parton It is interesting from the point of view of parametrizing the 



momentum flow out of a hadron vertex in a ^-independent way. However, the physics of a 
double distribution is equivalent to that of the corresponding off-forward distribution. 
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III. OFPD'S AND THE SPIN STRUCTURE OF THE NUCLEON 



One of the principal motivations for studying off-forward distributions is that they pro- 
vide a new class of observables for the internal structure of the nucleon. We believe that 
the physics of strong interactions is described by the fundamental theory-quantum chro- 
modynamics (QCD). However, QCD in the low-energy region is strongly interacting and 
highly relativistic, involving an infinite number of degrees of freedom. No theory of this 
type has ever been solved before with a clear physical insight. The only available tool now 
for solving QCD bound states is lattice QCD [|]. At this stage, the approach has not yet 
been developed to a degree of sophistication so that one can reliably calculate fundamental 
observables with respectable precisions. Therefore, any experimental information on the 
structure of the nucleon is valuable in helping to understand how nature constructed this 
particle of extreme importance to our existence. 

Since Gell-Mann and Zweig's quark model, many nucleon models have been proposed to 
correlate the experimental data from low and high energy probes. Although these models 
are not precisely QCD, they have been quite successful phenomenologically. However, the 
ultimate question confronting models is how the degrees of freedom used in them are related 
to those in QCD. This question is of high importance because experimental probes do not 
couple to model degrees of freedom. In practice, simple assumptions have been used in 
comparing model predictions and data, e.g., the constituent quarks are identified as QCD 
quarks at low-energy scales although a precise notation of the latter is unclear. These 
assumptions have been seriously challenged in 1987 when the European Muon Collaboration 
(EMC) measured, with an unprecedented precision, the proton's spin-dependent structure 
function Gi(Q 2 ,u) in polarized deep- inelastic scattering. Combining their data with the 
hyperon beta decay rates, augmented with the assumption of the flavor SU(3) symmetry, 
EMC extracted the fraction of the nucleon spin carried in the spin of quarks: 



This result is in flat contradiction with the cherished quark model prediction AS = 1. 
This so-called "spin crisis" has touched off an explosive activity in both theoretical and 
experimental communities. After ten more years of careful experimental studies, the original 
EMC result has essentially been confirmed, although the small x contribution to AS remains 
largely uncertain. A recent global analysis of the data can be found, for instance, in Ref. 



Many theoretical "solutions" to the "spin crisis" have been proposed in the literature. 
Apart from doubts about the experimental data and the interpretation of them, there have 
been serious attempts in explaining the EMC result from the fundamental theory. For 
instance, lattice QCD calculations have produced numbers close to the EMC result |33H35|, 
although the systematic effects of quenched approximation, finite lattice sizes, and large 
quark masses remain uncertain. Despite these, more fundamental approaches have not yet 
offered a complete picture about the spin structure of the nucleon. As to the naive quark 
model prediction, the key question is what has been alluded to abovel: namely, whether the 
quark model result can be directly compared with the experimental data, or whether the 
constituent quarks are the same as the QCD quarks probed in high-energy scattering. No 
convincing argument has been offered towards a positive answer. 




(29) 
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More recently, progress has been made in attempting to understand how QCD builds 
up the nucleon spin. The progress potentially has an important impact on other theoretical 
studies, and ultimately on our basic knowledge of the internal structure of the nucleon. On 
the one hand, it suggests a thorough examination of the spin structure in a fundamental 
way. On the other hand, they called for understanding the role of the gluon degrees of 
freedom in model building. As it turns out, a key element in unraveling the nucleon spin is 
the off-forward parton distributions. 

The history began with identifying the sources of the QCD angular momentum. Jaffe 
and Manohar wrote down the complete angular momentum operator in a free-field-theory 
form which contains naturally the quark and gluon spin and orbital contributions |!6| . The 
importance of the quark orbital angular momentum in QCD has been recognized before by 
Sehgal [B7|] and Ratcliffe in different contexts. The EMC result indicates, among others, 
that a large fraction of the nucleon spin is carried by other sources of angular momentum. In 
order to associate these contributions to possible observables, the QCD angular momentum 



operator has been rewritten in a gauge- invariant form [21 



^qcd — J q + J g 



(30) 



where 



Jo 



d cc x x Tfj 



Jn 



d 3 x 
•3 



d x x x (E x B) . 



(31) 



The quark and gluon parts of the angular momentum are generated from the quark and 
gluon momentum densities T q and E x B, respectively. E is the Dirac spin-matrix and the 
corresponding term is clearly the quark spin contribution. D = d + igA is the covariant 
derivative and the associated term is the gauge-invariant quark orbital contribution. 

With the above expression, one can easily construct a sum rule for the spin of the nucleon. 
Consider a nucleon moving in the ^-direction, and polarized in the helicity eigenstate A = 
1/2. The total helicity can be evaluated as an expectation value of J z in the nucleon state 



2 2 vw 



L q (n) + J q (fj) 



(32) 



where the three terms denote the matrix elements of three parts of the angular momentum 
operator in Eq. fl3"T|). The physical significance of each term is obvious, modulo the scale and 
scheme dependence indicated by fi. The scale dependence in £(/z) is generated from the U(l) 
axial anomaly, and there have been attempts to remove the scale dependence by subtracting 
a gluon contribution |39| , |40"f . Unfortunately, such a subtraction is by no means unique. Here 
we adopt the standard definition of AS(/i) as the matrix element of the multiplicatively 
renormalized quark spin operator. Note that the individual term in the above equation is 
independent of the momentum of the nucleon. In particular, it applies when the nucleon is 
travelling with the speed of light (the infinite momentum frame) [[y]] . 
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The scale dependence of the quark and gluon contributions can be calculated in pertur- 



bative QCD. By studying renormalization of the nonlocal operators, one can show f42| , |2T 





(33) 



As \i — > oo, there exists a fixed point solution 



Jq(oo) 



Jg((X)) 



3n 



f 



216 + 3n/ 
1 16 
216 + 3n/ 



(34) 



Thus as the nucleon is probed at an infinitely small distance scale, approximately one-half 
of the spin is carried by gluons. A similar result has been obtained by Gross and Wilczek in 
1974 for the quark and gluon contributions to the momentum of the nucleon |j4"3"| . Strictly 
speaking, these results reveal little about the nonperturbative structure of bound states. 
However, experimentally it was found that about half of the nucleon momentum is carried 
by gluons already at quite low energy scales [|7| . Thus the gluon degrees of freedom not only 
play a key role in perturbative QCD, but also are a major component of nonperturbative 
states, as they should be. An interesting question is then, how much of the nucleon spin is 
carried by the gluons at low energy scales? A solid answer from the fundamental theory is 
not yet available. Recently, Balitsky and this author have made an estimate using the QCD 
sum rule approach j44|. We found 



J 9 (H ~ lGeV) 



> e < uaGu >< uu > 



9 



(35) 



N 



which yields approximately 0.25. If this calculation indicates anything of the truth, the spin 
structure of the nucleon roughly looks like 



1 1 

- = 0.10(from -AS) + 0.15(from L q ) + 0.25(from J g ) 



(36) 



A quark model calculation of J g has also been made recently by Barone et al. ||45|| . The 
result also indicates J g ~ 0.25 at low-energy scales. 
By examining carefully the definition of J q ^ g 



p 1 - 

2 



d 3 x(x x T, 



Q,9J 



P 1 - 

2 



(37) 



one realizes that they can be extracted from the form factors of the quark and gluon parts 
of the QCD energy-momentum tensor T^ v . Specializing Eq. (Q) to (n = 2) 

(P'\T%\P) = U(P')[A q>g (t)^P . 

(38) 

Taking the forward limit in the n = component and integrating over 3-space, one finds 
that A qg (0) give the momentum fractions of the nucleon carried by quarks and gluons 
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(A q (0) + A g (0) = 1). On the other hand, substituting the above into the nucleon matrix 
element of Eq. (j37|) , one finds [21 



J q , 9 = ^ q ,M + B q ,M\ ■ (39) 

Therefore, the matrix elements of the energy-momentum tensor provide the fractions of the 
nucleon spin carried by quarks and gluons. There is an analogy for this. If one knows 
the Dirac and Pauli form factors of the electromagnetic current, Fi(Q 2 ) and i^Q 2 ), the 
magnetic moment of the nucleon, defined as the matrix element of (1/2)/ d 3 x(x x j) z , is 
Fi(0) + F 2 (0). 

Since the quark and gluon energy-momentum tensors are just the twist-two, spin-two, 
parton helicity-independent operators, we immediately have the following sum rule from the 
off-forward distributions 

J dxx[H g (x, t) + E g (x, f , *)] = A q (t) + B q (t) , (40) 

where the £ dependence, or C q {t) contamination, drops out. Extrapolating the sum rule to 
t = 0, the total quark (and hence quark orbital) contribution to the nucleon spin is obtained. 
A similar sum rule exists for gluons. 

Thus a deep understanding of the spin structure of the nucleon may be achieved by 
measuring the OFPD's from high-energy experiments. 

IV. WHAT DO WE KNOW ABOUT OFPD'S? 

Comparing with the ordinary parton distributions, off-forward distributions depend on 
two additional variables: the t-channel longitudinal momentum fraction £ and the invariant 
mass t. Because of these, the OFPD's contain much more information and are necessarily 
more difficult to model in practice. At present, there is little experimental data to constrain 
these distributions. In this Section, we will describe a few of the recent studies of the 
distributions at low-energy scales where they reflect directly the nonperturbative structure 
of the nucleon. We will also discuss their evolution equations which have been computed to 
leading logarithmic order in perturbation theory. Some qualitative features of the evolution 
to an asymptotically large energy scale will be highlighted. 

A. OFPD's at Hadron-Mass Scale 

A few rigorous results about the OFPD's are known. First of all, in the limit of £ — > 
and t — > 0, they reduce to the ordinary parton distributions. For instance, 

H q (x,0,0) = q(x) , 

H q (x,0,0) = Aq(x) , (41) 

where q(x) and Aq(x) are the unpolarized and polarized quark densities. Similar equations 
hold for gluon distributions. For practical purposes, in the kinematic region where 
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t\ « M N and £ « x (42) 

an off-forward distribution may be approximated by the corresponding forward one. The first 
condition, <J\b\ <<C M N , is crucial - otherwise there is a significant form-factor suppression 
which cannot be neglected at any x and £. For a given t, £ is restricted to 



K| < v/-t/(M2 - t/4) . (43) 

Therefore, when J\t\ is small, £ is automatically limited and there is in fact a large region 
of x where the forward approximation holds. 

The first moments of the off-forward distributions are constrained by the form factors of 



the electromagnetic and axial currents. Indeed, by integrating over x, we have |21 



f^dxH q {x,^t)=F q l {t) 
J dxE q (x,S,t)=Fg(t) 
j\xH q {x,£,t)=G q A (t) 



i 



dxE q (x,^t) = G q P (t) . (44) 



The t dependence of the form factors are characterized by hadron mass scales. Therefore, it 
is reasonable to speculate that similar mass scales control the t dependence of the off-forward 
distributions. In the small x region, it has been argued that the t dependence corresponds 
to a mass scale 4m 2 or something close to the inverse slope of the elastic PP cross section, 
1/B = 0.1 GeV 2 g§. 

Martin and Ryskin made a very important observation that one can obtain an 
upper bound on the off-forward distributions by using the Schwarz inequality [f46 
For instance, consider the quark distribution in Eq. ( O) and take a ket 



bx ((x + 0P + , k ± ) \P)/^x~Tl - 6 A ((x - 0P + , k± + A ± ) |P')/V^£. Since > 0, we 

have the following constraint 

F q (x,0<l(q(x 1 ) + q(x 2 )) , (45) 

for x > £, where x\ = (x + £)/(l + £) and x 2 = (x — £)/(l — £) which is different from that 
in Eq. (pBl). Similarly, one has 

- F q (x, < i (q(- Xl ) + q(-x 2 )) , (46) 

for x < — £. In the middle region, however, the vacuum subtraction ruins a useful constraint. 
An upper bound for the off-forward gluon distribution can be derived in a similar way: 

2xF g (x, £) < - (xisf(xi) + x 2 g{x 2 )) , (47) 

for a; > £. Because of the different definitions of x 2 , the above inequality is slightly different 
from Martin and Ryskin's. However, they are practically the same in the small x region. It 



16 



must be cautioned that the above inequalities are derived by disregarding renormalization 
of the distributions, and hence are correct in the leading-logarithmic approximation only. 
A first model calculation of the quark off-forward distributions was reported in Ref. |47| . 



The calculation was done in the MIT bag model, whose parameters are adjusted so that 
the electromagnetic form factor and the Feynman parton distributions are best produced. 
The shape of the distributions as a function of x are rather similar at different t and £: 
peaked around x = 0.2 ~ 0.4, and tapering off rather quickly as x — > 1, and negative, 
where the calculation cannot be trusted. The t dependence of the energy-momentum form 
factors is controlled by a mass parameter between 0.5 and 1 GeV 2 . The £ dependence of the 
distributions is surprisingly weak; a clear physical reason is unknown. 

The off-forward quark distributions were also studied in the chiral quark-soliton model by 
Petrov et al. |48|]. In contrast to the bag model results, the chiral soliton model yields a rather 



strong £ dependence of the off-forward distribution. The model also predicts qualitatively 
different behaviors in the regions \x\ > £ and |x| < £, with discontinuities at x = ±£. 
Although the discontinuities are believed to be an artifact of the calculation, the qualitatively 
different features in different regions are in line with the interpretation of the OFPD's. 
From Eq. (|T4j), we understand that the off-forward distributions have similar physics as the 
forward distributions in \x\ > £. When \x\ < £, however, the distributions reflect t-channel 
exchanges of meson and glueball quantum numbers. Therefore, a physics-oriented modelling 
in this region may consider the convolution of the meson cloud of a nucleon with the leading 
light-cone wave function of the mesons. 

Models cannot be used yet to calculate off-forward gluon distributions in the small x 
region, which are phenomenologically important 0-01 . It is believed in the literature that 
the distributions in this region and negligible t are independent of £ [22,14]. As we mentioned 
earlier, this might be correct as long as £ <C x. For £ ~ x, however, the best argument so 
far is offered by Martin and Ryskin ||46|| . From the diffractive scattering phenomenology, 
they argue that the distribution F g (x,£) at small x should be bounded below by the forward 
distribution at x 2 = {x + 0/(1 + 5 

2xF g (x,0>x 2 g(x 2 ) . (48) 

Combining this with the upper bound in Eq. p7|), they concluded phenomenologically when 
x < 10" 2 , 

2xF g (x, f , fi) ~ x 2 g(x 2 , fi) , (49) 

at /j 2 = 0.4 GeV 2 for the GRV gluon distribution § and at /j 2 = 1.3 GeV 2 for the MRS(R2) 
distribution 101. 



B. Evolution Equations 

Parton distributions are renormalization scale and scheme dependent because the defin- 
ing operators are. Thus, a natural starting point for discussing scale evolution is to consider 
renormalization of twist-two operators. However, it is well known that the parton evolu- 
tion has a simple physical interpretation in light-cone coordinates and gauge, as exemplified 
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by the DGLAP evolution kernels. Therefore, we shall also consider the off-forward parton 
evolution kernels. 

Renormalization of the twist-two operators appearing in Feynman parton distributions 
is a special example of a more general case. For a given twist-two operator of spin n shown 
in Eq. ([I]), apart from its mixing with the gluon operator of same quantum numbers it also 
mixes with the following operators with total derivatives, 



— <->M2i+l 



id^ 1 ■■■id m ij;iV 



■iV 



(50) 



Consideration of such mixing is already necessary in studying the evolution of the leading- 
twist meson wave functions. The answer was first obtained by Efremov and Radyushkin 



17], and by Brodsky and Lepage [27]. Actually, at the leading logarithmic order the answer 



may be guessed from the naive conformal symmetry which is broken by quantum effects only 



at the next-to-leading-logarithmic order |50|J5T[ . The following combination of the twist-two 
operators furnishes a representation of the special conformal symmetry group, 



6 n = {id ■ ny-^c^ 



i T> -n 

id+i~d. 



(51) 



where the C^^ix) are the Gegenbauer polynomials of order 3/2. Its evolution takes exactly 



the diagonal form of that for the operator without the total derivatives [[17 



where 



7n = 2C F 



n 1 

4 ^7- 3 - 

. i=l b 



"In 

2/3 ~ 



O n (fi 2 



(52) 



n(n + 1 



(53) 



and Po = 11 — 2iif/3. Cp = 4/3 for the SU(3) color group and rif is the number of light 
quark flavors. For the leading-twist pion wave function, the above result leads to a general 
expansion in terms of Gegenbauer polynomials. 

Generalization to the singlet case is in principle straightforward. For the twist-two 
operators with natural parity, the evolution was first worked out by Chase |52][ . For those 
with unnatural parity, it was worked out by Ohrndorf [53|, Shifman and Vysotsky j53 



and Baier and Grozin f5"5(| . The conclusion of these studies is that the gluon towers of 
operators evolve in the leading-logarithmic approximation exactly like the forward case if 
they are combined according to the Gegenbauer polynomials of order 5/2. These results 
were confirmed later in a more elegant form of light-ray or string operators by Geyer et al. 
| 56| , and by Balitsky and Braun [|57 . 



In light of these developments, evolution of off-forward distributions can simply be worked 
out by converting the anomalous dimensions of the composite operators into evolution ker- 



nels. However, as in the case of the DGLAP evolution fllql , one can obtain the same results 
more directly through studying parton splitting processes. A first such study for the non- 
forward, helicity-independent parton evolution was made by Bartels and Loewe in the small 
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x limit ||. 
glet mixing 



Gribov et al. presented the complete leading-order kernels including the sin- 
ToU . Miiller et al. found the non-singlet evolution kernel by generalizing the 
More recently, the issue has been reexamined from different 



Brodsky-Lepage result JT9 



angles by Radyushkin |23|,[3(J], Ji |58|, Chen |59 ], Frankfurt et. al. [6(J, Bliimlein et al. |6l 



As an example of the off-forward evolution, we cite below the 



and Martin and Ryskin PB] 
helicity-independent result in Ref. 

According to Eq. (TI^), it is natural that the distributions in three different kinematic re- 
gions evolve differently. In the region x > £ where we have a quark creation and annihilation, 
the evolution equation is, 



D Q F NS (x,£,Q 2 ) a s {Q 2 



where 



D ( 



DlnQ 2 



d 



2tt 



dy 



P> 



NS 




g s {Q 2 

D\nQ 2 d\nQ 2 2n 
The parton splitting function is 



-c. 



3 

2 + J, u 



dy 



x 



dy 



-e y 



X 



te 



P NS (x,£,e) = C F - 



x 2 + 1 - 2e 



1 -x + ie)(l -£ 2 ) 



(54) 



(55) 



(56) 



The end-point singularity is cancelled by the divergent integrals in Dq/ DhiQ 2 . Obviously, 
when £ = 0, the splitting function becomes the usual DGLAP evolution kernel ||16[ . For 
— £ < x < £ where we have creation or annihilation of a quark- ant iquark pair, the evolution 
takes the form, 



D Q F NS (x,£,Q 2 ) _ a s {Q 2 



DlnQ 2 



2tt 



1 dy , ( x £ e 

*NS \ ~ ' _ ' _ 

* y \v v v, 



dy pl fx £ e 

*NS 1—1 _ ' _ 



-i y 



y y y, 



where 



Pns( x i £> e ) 



2C F 



x + £ 



2? 



x + it 



F NS (y,Z,Q 2 ) , 

(57) 



(58) 



When £ = 1, the kernel reduces to the Brodsky and Lepage kernel P? |. For x < — £ where 
we have antiquark creation and annihilation, the evolution takes the same form as Eq. 
apart from the replacement J^ 1 — ► — /f v 

The helicity-dependent off-forward distributions were introduced in Ref. |58|] and their 
evolution was also studied in that paper. Subsequently, the kernels have been confirmed 
by the studies of Bliimlein et al. |)T and Balitsky and Radyushkin [JS^] . The helicity-flip 
off-forward quark distributions were introduced by Collins et al. f25| , and their evolution was 

The helicity-flip gluon distributions were defined by 
who also derived the leading-logarithmic evolution. It is worthwhile 



studied by Belitsky and Miiller |63[ 
Hoodbhoy and Ji ]2?_ 

to point out that the combinations of the helicity-flip quark and gluon operators according 
to the Gegenbauer polynomials of order 3/2 and 5/2 respectively also form representations 
of the naive conformal symmetry group and hence have diagonal evolutions at the leading- 
logarithmic level. 
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C. Effects of Evolution 



As in the forward case, the off-forward parton evolution equations can be solved in 
many different ways. A few approaches have already been suggested in the literature 
30| , |60|j6^ , |65|^6| . One important aspect of the evolution is the asymptotic form of the 

(also 



off-forward distributions. Radyushkin found 
distributions, 



55|) for the non-singlet quark 



F NS (x,£,t,Q 2 



? (1 



e 



(59) 



where F(t) is the appropriate electromagnetic form factor. This asymptotic form has inter- 
esting implications. 

When £ = 0, the asymptotic form is a 5-function in x. This reminds us of the well-known 
fact that as Q 2 —>■ oo, parton distributions migrate to smaller and smaller x. For £ = 1, 
the evolution is that of a meson wave function with an asymptotic limit 1 — x 2 , which has a 
broad shape in x. Therefore, it is quite plausible that for an intermediate £, one would see 
the migration of the distribution from large to small |x| if \x\ > £. However, in the middle 
region |x| < £, the total strength increases and at the same time, the distribution evens out 
toward 1 — x 2 . These features of evolution were observed in the numerical calculations by 
Belitsky et al. |)4] and Mankiewicz et al. [65 



The evolution effects for the off-forward gluon distribution at small x are especially 
interesting from the phenomenological point of view. These effects have been studied by 



Frankfurt et al. |60], and Martin and Ryskin [46]. In Ref. |]60 |, they assumed negligible off- 
diagonal effects at Q = 1.6 GeV, and then evolved the off-diagonal distribution to Q — 7, 
17, 41, and 110 GeV at x x = 1.1 x 10~ 2 , 1.1 x 10~ 3 , and 1.1 x 10" 4 . For x 2 /x x = 10~ 2 
to 1, they found that the ratio of off-diagonal to diagonal distributions ranges from 1.2 to 
1.7. They pointed out that for a fixed x 2 — Xi, the off-diagonal distribution at large Q is 
obtained from that in the much larger x region, however, with the same difference x 2 — x\ 
where the off-diagonal effect is very small. Therefore at very small x, unless Q is very close 
to hadron mass scales, the off-diagonal effects come mostly from the perturbative evolution. 
These effects are less than a factor of 2 for the kinematics region considered. 

Similar studies have also been made in Ref. [[46]. The qualitative conclusion is the 



same although quantitatively there are differences. Martin and Ryskin pointed out that the 
off-diagonal evolution happens faster than the diagonal one because the difference in the 
splitting functions is positive. Furthermore, the difference is only important when Q is close 
to the final target scale and x becomes comparable to the size of £. 



V. PROBING OFPD'S IN HIGH-ENERGY EXPERIMENT 



Eventually, usefulness of the OFPD's depends on whether they can actually be measured 
in any experiment. The simplest, and possibly the most promising, type of experiments are 
hard electro- or deep-inelastic production of photon, mesons, and perhaps even lepton pairs. 
In this section, we will consider separately two experiments: deeply virtual Compton scat- 
tering (DVCS) in which a real photon is produced, and diffractive meson production. There 
are practical advantages and disadvantages from both processes. Real photon production 
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is, in a sense, cleaner but the cross section is reduced by an additional power of a em . More- 
over, the DVCS amplitude interferes with the Bethe-Heitler one. Meson production may be 
easier to detect, however, it has a twist suppression, l/Q 2 . In addition, the theoretical cross 
section depends on the unknown light-cone meson wave function. At present, there exist 
already some data from HERA on meson production flBUfl . More experiments will be done 
in the future at HERA, CERN, and Jefferson Lab. 



A. Deeply Virtual Compton Scattering 



Deeply virtual Compton scattering was proposed in Ref. [2T| as a practical way to mea- 
sure the off-forward distributions. The result was confirmed in Ref. |yj using the double 
distributions. The general two-photon processes had actually been studied earlier by Watan- 
abe || and Miiller et al. | l9f , and later by Chen [59]. Recently, there have been a number 
of new results obtained about this process, which will be described below. 

Consider the virtual photon scattering shown in Fig. 2a, where the momenta of the 
incoming (outgoing) photon and nucleon are q(q') and P(P'), respectively. The Compton 
amplitude is defined as 



d^ze^ ( P' 



tj^ — \r 



p 



(60) 



27 " V2, 

where q = (q + q')/2. In the Bjorken limit, — q 2 and P ■ q —>■ oo and their ratio remains 
finite, the scattering is dominated by the single quark process shown in Fig. 2b in which a 
quark absorbs the virtual photon, immediately radiates a real one, and falls back to form 
the recoiled nucleon. In the process, the initial and final photon helicities remain the same. 



The leading-order Compton amplitude is pl|j 




J2e 2 q F q (x,^t,Q 2 



x + £ — ie 



Y,e 2 q F q (x^^Q 2 



(61) 



where n and p are the conjugate light-cone vectors defined according to the collinear direction 
of q and P, and g^ 1 is the metric tensor in the transverse space. £ is related to the Bjorken 
variable x B = -q 2 /(2P ■ q) by x B = 2£/(l + f )• 





(a) (b) 
FIG. 2. a). Compton scattering; b) Leading Feynman diagrams for DVCS. 
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The same DVCS final state can also be produced through the Bethe-Heitler process in 
which the initial or final state electron (or muon) radiates a real photon and at the same 
time scatters elastically off the target nucleon. This process can be calculated as accurately 
as the data on the nucleon elastic form factors. Physically, the Bethe-Heitler process is a 
background to DVCS and can overshadow the latter signal if the former is too large. On the 
other hand, with an appropriate size of Bethe-Heitler amplitude, one may hope to measure 
its interference with the DVCS amplitude which then can be directly extracted from data. 
The complete leading-order DVCS cross sections and the interference with the Bethe-Heitler 
amplitude can be found in Ref. 

Recently, Guichon made some interesting estimates of the cross sections at COMPASS 
and Jefferson Lab energies |67]. He found that while at low scattering energy, the Bethe- 



Heitler background is quite large, it drops significantly at high energies when the virtual 
photon flux is large. In fact, if one is focused on production of real photons along the 
direction of the virtual photon, the DVCS cross section is quite significant at COMPASS 
energy, and certainly far above the Bethe-Heitler background. As far as studying the orbital 
angular momentum of quarks is concerned, this is the most interesting kinematic region. 
Guichon's estimate is encouraging for a real experiment. 

The one-loop corrections to DVCS have recently been studied by Osborne and Ji pS 



They have also been studied by Miiller |69|] and by Belitsky and Miiller [ 33 1 using the 



constraints from conformal symmetry. Results of these studies have been confirmed by 



Mankiewicz et al. |70|. Besides their obvious use in precision analysis, the results indicate 
that DVCS is factorizable at the next-to-leading order. A higher-order calculation of the 
coefficient functions in the renormalon-chain approximation has been done by Belitsky and 



Schafer JTl]. An all order proof of the DVCS factorization was first given by Radyushkin 



and recently in a different perspective by Osborne and Ji |]72|| , and Collins and Freund 
The simple DVCS scattering mechanism can be tested through scaling relations and 



selection rules, as in the case of deep-inelastic scattering. In Ref. ||74|| , Diehl et al. considered 
the possibility of using <p dependence of the cross section to test the single quark scattering 
picture, where </> is the angle between the lepton and hadron scattering planes. They observed 
that through QCD radiative corrections, the photon helicity-flip amplitude may contribute 
to DVCS. A recent study by Hoodbhoy and Ji showed that this indeed happens |p8|| . The 
key element here is the helicity-flip off-forward gluon distributions, which decouple in the 
forward limit due to angular momentum conservation. They pointed out that a nonvanishing 
leading-twist photon helicity-flip amplitude signals unequivocally the presence of gluons in 
the nucleon. 

There are other processes that are as clean as the DVCS, but more difficult to mea- 
sure experimentally. For instance, Z° production was considered in Ref. |§. Since Z°'s are 
signaled through lepton pairs, one can also consider lepton-pair production through inter- 
mediate time-like virtual photons. Theoretical treatment of these is exactly the same as for 
DVCS. 
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B. Vector Meson Production 



Heavy quarkonium production was first considered by Ryskin as a way to measure the 
Feynman gluon distributions at small x [[□]]• A leading-order diagram for the process is 
shown in Fig 3a, where the virtual photon fluctuates into a cc pair which subsequently 
scatters off the nucleon target through two-gluon exchanges. In the process, the pair transfers 
away certain amount of its longitudinal momentum and reduces its invariant mass to that 
of a J/ty. The cross section is, 



da 
~dt 



(7* + P -> J/tf + P') 



i67r 3 Mr e - 

3a cm Q e 



2\1 2 



(62) 



where Q 2 = (Q 2 + M 2 )/4, M is the J/ip mass, and T e + e - is the decay width into the lepton 
pair. The equation was derived in the kinematic limit s ^> Q 2 ^> M 2 3> t and the Fermi 
motion of the quarks in the meson was neglected. Two other important approximations 
were used in the derivation. First, the contribution from the real part of the amplitude 
is neglected, which may be justifiable at small x. Second, the off-forward distributions 
are identified with the forward ones. As we have discussed in IV.A, even for small t, the 
approximation is good only if the off-diagonal evolution effects can be neglected. A more 



refined treatment of the process and comparison with HERA data can be found in Ref. |T3 





FIG. 3. a). A representative Feynman diagram for neutral meson production through 
off-forward gluon distributions in the nucleon; b) A representative Feynman diagram for meson 
production through off-forward quark distributions in the nucleon. 

The above result was extended to the case of light vector-meson production by Brodsky 
et al., who considered the effects of meson structure in perturbative QCD [|12j. They found 
a similar cross section, 

^ N ^ VN ) = ^-qs > (63) 



da 



t=o 



where the dependence on the meson structure is in the parameter 

1 f dz aV( ' ' ~ 

2 J {z 



Vv 



dzcf) 



v 1 



(64) 



and (f) v (z) is the leading-twist light-cone wave function. Evidently, the above formula reduces 
to Ryskin's in the heavy-quark limit (<f) v (x) = 5(x — 1/2)). More discussions on vector 
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meson production along this line, including comparison with recent data, can be found in 
Ref. 



13,14,75 



Radyushkin was the first to calculate the amplitude of hard diffractive electroproduction 
in terms of off- forward gluon distributions |23|j3(| . With the virtual photon and vector meson 
both polarized longitudinally, his result corresponds to 



da 
~dt 



( 7 *iV -> VN) 



A7rT v m v a 2 s (Q)rjl 



2x 



B 



dx 



1 



IE 



X 



F g (x,£,t) 



(65) 



where again xb = 2£/(l + £). The above formula is valid for any xb and t smaller or around 
hadron mass scales. Hoodbhoy has also studied the effects of the off-forward distributions 
in the case of J/ip production p4| . He found that Ryskin's result needs to be modified in a 
similar way once the off-forward effects becomes important. 

Collins, Frankfurt, and Strikman have made an extensive analysis of the high-order 
QCD contributions to hard diffractive meson production |25[| . Using the standard tools 
described in 0, they showed that the leading- twist contributions are factorizable to all 
orders in perturbation theory. The factorization theorem provides a rigorous QCD basis for 
calculating the process in terms of the off- forward distributions. In fact, the paper outlined 
a recipe for calculating a whole class of processes in perturbation theory. Mankiewicz, Piller 
and Weigl recently calculated the neutral pseudoscalar and vector meson production cross 
sections in terms of the non- forward quark and gluon distributions |70| . Shown in Fig. 3b is 
a Feynman diagram in which meson production depends on a quark distribution. They also 
investigated charged particle production in which the initial nucleon has a different charge 
from the final one JTIJ . The charge- changing off- forward distributions can be obtained from 
the charge- conserving ones using isospin symmetry. The result can be generalized to the 
octet baryons if the SU(3) flavor symmetry is good. One can also consider production 
of heavy-flavored baryons from the nucleon target. Unlike the heavy-quark distributions in 
light hadrons, the corresponding off-forward distributions are not suppressed by heavy-quark 



masses. 



VI. SUMMARY AND OUTLOOK 

The off-forward parton distributions are a new type of nucleon observables which have 
recently generated much interest in hadron physics community. On the one hand, the 
OFPD's contain all the form factors of the twist-two operators which generalize the well- 
known electroweak currents. The form factors of the spin-two, twist-two operators (energy- 
momentum tensor) are especially interesting in unravelling the spin structure of the nucleon. 
On the other hand, they summarize succinctly the leading structural dependence in a class 
of high-energy processes involving the nucleon. The partonic interpretation of the OFPD's 
is quite simple- it represents parton-nucleon scattering amplitudes on the light-cone. 

Recent studies of the OFPD's can largely be classified into three categories: modelling 
the new distributions, studying their scale dependences, and investigating processes that 
are sensitive to them. In terms of model studies, our accomplishments so far are rather 
limited. Even though there are bag and chiral soliton model calculations, these calculations 
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need much improvement in the future. In particular, the physics in the \x\ < £ region has 
not been entirely captured in these calculations. The models are yet to be extended to 
including gluon degrees of freedom which are known to be important even at low energy 
scales. Indeed, a vanishing helicity-flip gluon distribution at low-energy scales will remain 
zero at other scales because it has no mixing with any quark distribution. A number of 
rigorous constraints discussed in Section IV.A can be helpful in modelling: Apart from the 
special limits and n — 1 moments, there exist useful upper bounds for the OFPD's. 

There have been many studies in perturbative evolution. At the leading-logarithmic level, 
evolution equations now have been obtained for all leading-twist OFPD's. At next-to-leading 
order, one needs two-loop anomalous dimensions which exist only for the spin-independent 
non-singlet operators |77| . Calculations of additional two-loop anomalous dimensions will be 
needed as data on the OFPD's become available. Recent works of Miiller and Belitsky illus- 
trate the power of exploiting the anomalous conformal symmetry to obtain these anomalous 
dimensions J7S|. Numerical studies of evolution effects have just begun. Qualitatively we 
can understand the evolution effects from its asymptotic behavior. Quantitatively, one needs 
to establish high-precision, high-speed evolution codes which allow efficient data fitting in 
the future. 

The leading-order cross sections for DVCS and meson production can be found in the 
literature. For DVCS, the cross section has also been computed to the next-to-leading order 
in a s . A similar study for meson production can be made, and will be necessary once 
precision data become available. Higher-twist effects are important at low Q 2 and must be 

investigated thoroughly. Generally speaking, they are more complex to study than that 
in deep-inelastic scattering because of the extra photon or meson, and the extra soft scale 
t. Additional processes may be explored to probe the OFPD's, including ones in which the 
final state baryon may contain strange or heavy flavors. Polarization observables need to 
be seriously investigated. A solid experimental confirmation of the scattering mechanisms 
involving the OFPD's is itself very significant. It involves testing scaling relations and 
selection rules. 

Finally, we notice that the OFPD's are functions of multi- variables. As such, it would 
be difficult to get a complete picture of them in a short time period and in a single type 
of experiments. Moreover, experimental cross sections cannot be used to extract the distri- 
butions directly; they usually involve integrals over parton momenta. Therefore, it will be 
necessary to parametrize the OFPD's and fit the parameters to data. In finding appropriate 
ways to make parametrizations and fits, model studies will play important role, unless of 
course lattice calculations can make a significant breakthrough in the near future. 

Note Added: After this paper was submitted for publication, the author has noticed 
some new works in the literature which are relevant to the subject of this paper. In Ref. 
f79fl , Radyushkin studied both hadron form factors and wide-angle Compton scattering using 
nonforward parton densities. In Ref. [j50|] , diffractive dijet photoproduction was proposed as 
a probe of the off-diagonal gluon distribution. In Ref. f81| , certain positivity constraints for 
off-forward parton distributions were derived. In a different line of development, the orbital 
angular momentum distribution in the nucleon was investigated |82| . 
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